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Abstract 

In the Randah-Sundrum cosmological models, a (3+l)-dimensional brane subject to a Z2 orb- 
ifold symmetry is embedded in a (4+l)-dimensional bulk spacetime empty except for a negative 
cosmological constant. The unperturbed braneworld cosmological solutions, subject to homo- 



' geneity and isotropy in the three transverse spatial dimensions, are most simply presented by 

■ means of a moving brane description. Owing to a generalization of Birkhoff 's theorem, as long 



as there are no perturbations violating the three-dimensional spatial homogeneity and isotropy. 



I the bulk spacetime remains stationary and trivial. For the spatially flat case, the bulk spacetime 

is described by one of three bulk solutions: a pure AdS^ solution, an AdS'^-Schwarzschild black 
hole solution, or an yldS'^-Schwarzschild naked singularity solution. The brane moves on the 



boundary of one of these simple bulk spacetimes, its trajectory determined by the evolution of 
the stress-energy localized on it. We derive here the form of the Israel matching conditions for 



■ the linearized cosmological perturbations in this moving brane picture. These Israel matching 

conditions must be satisfied in any gauge. However, they are not sufficient to determine how to 
describe in a specific gauge the refiection of the bulk gravitational waves off the brane boundary. 
In this paper we adopt a fully covariant Lorentz gauge condition in the bulk and find the sup- 
plementary gauge conditions that must be imposed on the boundary to ensure that the reflected 
waves do not violate the Lorentz gauge condition. Compared to the form obtained from Gaussian 
normal coordinates, the form of the Israel matching conditions obtained here is more complex. 
However, the propagation of the bulk gravitons is simpler because the coordinates used for the 
background exploit fully the symmetry of the bulk background solution. 



1 



I. INTRODUCTION 



This paper discusses gravitational gauge invariance and gauge fixing for tlie complete def- 
inition of the problem of computing the cosmological perturbations in a one-brane Randall- 
Sundrum scenario . In this cosmological scenario, the (4+l)-dimensional bulk space- 

time consists of a semi- infinite region, empty except for a negative cosmological constant, 
bounded by a (3+l)-dimensional timelike surface on which a singular Z2-symmetric distri- 
bution of stress-energy is localized. The timelike boundary moves in the semi-infinite bulk 
spacetime, its trajectory determined by the evolution of the stress-energy on it. We assume 
that the bulk spacetime metric is obtained from a small linearized deviation from an exact 
solution to the Einstein equations, so that Qat = g^ab + ^ab in the bulk. 

In this paper we establish the form of the boundary conditions coupling the degrees of 
freedom on the brane to those in the bulk when Lorentz gauge is chosen in the bulk. For 
the simplest situations such as a Minkowski or de Sitter brane, this is not the simplest 
gauge choice. For these special situations it is advantageous to employ Gaussian normal 
coordinates, so that h^^ and h^^ vanish. However, for considering realistic cosmological 
backgrounds, where the universe on the brane has an essentially arbitrary expansion history, 
Gaussian normal coordinates are particularly awkward, and in many cases develop artificial 
coordinate singularities. For studying more general, realistic cosmological solutions, it is 
advantageous to employ instead a "moving brane" description, where the unperturbed bulk 
remains static and trivial, described by coordinates such as the Poincare coordinates 



and the brane "moves," tracing out a timelike trajectory described as z = Zb{t). The mov- 
ing brane description has been developed in |^^J, in contrast to the fixed brane approach 
used in To consider perturbations in this moving boundary description, it is advanta- 
geous to adopt a "covariant" gauge in the bulk, exploiting the full AdS^ symmetry of the 
unperturbed bulk. One such gauge is Lorentz gauge, where 



In this gauge h^^ and /i^s no longer necessarily vanish and the form of the linearized match- 
ing condition contains additional terms. Moreover, five additional auxilliary ("gauge") 
boundary conditions must be imposed at the brane. One might say that these are the 



ds'^ = — dz^ — dt^ + dx\ + dx\ + dx\ 



(1) 




(2) 
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boundary conditions for the five longitudinal ("pure gauge") modes present in (4 + 1)- 
dimensional Lorentz gauge. The auxilliary gauge boundary conditions must be chosen in 
such a way that no forbidden polarization components are emitted or reflected into the 
bulk from the brane boundary. 

In this paper we establish the form of the linearized Israel matching condition and 



auxilliary gauge conditions for Lorentz gauge. In a forthcoming paper p3| we apply the 
formal results developed here to the actual problem of computing the evolution of the 
cosmological perturbations of the coupled brane-bulk system for various expansion histories 
and various choices for the physical degrees of freedom on the brane. 

Linearized braneworld cosmological perturbations have been previously studied by a 
large number of authors. (See refs. [t^ jssl for a partial sampling of the literature.) The 
majority of this work relies on the pure scalar-vector-tensor decomposition in the three 
transverse spatial dimensions, exploiting the three-dimensional homogeneity and isotropy 
of the background solution, so that the problem problem may be decomposed into separate 
sectors that are either pure scalar, pure vector, or pure tensor with respect to the three 
transverse spatial dimensions. As explained in Q and 3, generalizing on previous work 
in 3], each such sector contains only a single harmonic with respect to the three trans- 
verse spatial dimensions and its dynamics are described by a (H-l)-dimensional partial 
differential equation in the yt-plane. A sector corresponding to a pure tensor harmonic is 
described by a single coefficient function defined on the yt-plane, of scalar character with 
respect to y and t, whose evolution is described by a (1 + l)-dimensional partial differential 
equation. A sector corresponding to a pure vector harmonic has coefficient functions of 
both scalar and vector character with respect to y and t; however, as Mukohyama 0] and 

can be expressed in terms of a 
single "master" scalar function of pure scalar character with respect to y and t. Similarly, 
for the sector corresponding to a pure scalar harmonic, the resulting coefficient functions 
are of scalar, vector, and tensor character with respect to y and t. However, these in turn 
may be re-expressed in terms of another single "master variable" of scalar character. The 
approach adopted in this paper differs in that the gauge is fixed in a way that is completely 
local and independent of any particular choice of background coordinates or foliation of the 
background spacetime. It should be noted that the pure scalar-vector-tensor decomposition 
presupposes a high degree of symmetry and is nonlocal. 

The organization of the paper is as follows. In section |n] we analyse some simplified 
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examples that can be solved by expanding into plane waves in the bulk and individually 
matching at the boundary. In particular, we consider reflection of electromagnetic waves off 
a plane perfectly conducting wall under Lorentz gauge, and also the reflection of linearized 
gravity waves off a planar "gravitational mirror" embedded in a ffat (i.e., Minkowski) 
background. A "gravitational mirror" is a Z2-symmetric brane having a vanishing stress- 
energy at zeroth order. Section IIIII derives the equation of motion of the bulk metric 
perturbations in a general background in Lorentz gauge and characterizes the residual gauge 
freedom and its relation to reparameterization invariance on the brane. Section HVl considers 
the generalization of the results of section |n] to branes curved at zeroth order embedded in 
bulk backgrounds curved at zeroth order. Section |V| discusses stress-energy conservation on 
the brane and derives its expression at linear order. In section IVII the general form of the 
perturbative Israel matching condition and its divergence is presented. This result is then 
applied to show that the auxilliary boundary conditions proposed in section IIVI are in fact 
consistent with Lorentz gauge. Finally, section IVIII summarizes the results in a concise and 
self-contained form. Details of the derivation of the evolution equations of the bulk metric 
perturbations and the consistency of Lorentz gauge in the bulk, of the linearized Israel 
matching condition, and of the extrinsic curvature perturbations have been, respectively, 
relegated to three appendices. 

Finally, we make explicit the notational conventions used in this paper. We use the 
metric signature (— , -|-, +, +). Greek indices /x, z/, . . . denote 4-vectors in (3+1) dimensions; 
lowercase latin indices denote 5- vectors in a (4+l)-dimensional bulk spacetime (and in 
some cases, obvious from the context, one of arbitrary dimension as well). The uppercase 
latin indices A,B,... denote directions parallel to a brane on the boundary of the bulk 
spacetime, and denotes the inward unit normal with respect to the boundary. The 
order of the derivatives with the semicolon convention is as in the example, Wa-bc = 
A'^C^B'^'VcVbWa- The semicolon ; shall denote the covariant derivative with respect to 
the bulk and | denotes the covariant derivative with respect to the brane. We adopt the 
convention W^-bc — W"'.cb = R"'dcbW'^ for the Riemann tensor. On the boundary, we shall 
frequently use Fermi normal coordinates with respect to the directions on the surface of the 
brane, continuing these coordinates into the bulk using the Gaussian normal prescription. 
In this way the connection coefficients vanish at a given point, but in general not their 
partial derivatives. We shall also define = (^"^V^) for the bulk covariant derivative and 
define to be the corresponding derivation where the tilde indicates that the connection 
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is with respect to the induced metric on the brane. We shall also sometimes suppress raising 
and lowering of indices, it being implied that this is accomplished using the appropriate 
zeroth order metric. Accordingly the Einstein convention of summing over repeated indices 
is implied even when the pairs of repeated indices are both upper or both lower for more 
orderly looking expressions. We set = (SttG) where G is Newton's gravitational constant. 

II. SOME FLAT SPACE EXAMPLES: AN ACCELERATING CONDUCTOR 
AND PLANE GRAVITATIONAL MIRROR IN MINKOWSKI SPACE 

In this section we consider some very simple illustrations of what is to follow in which 
electromagnetic and linearized gravitational waves reflect from a flat planar boundary em- 
bedded in flat Minkowski space. In these simple cases the appropriate boundary conditions 
can be established merely by considering a decomposition into plane waves in the bulk. 
Later in the paper, when curved boundaries and a curved bulk are considered, the tech- 
nique of expanding into plane waves can no longer be exploited. Rather it is necessary to 
modify slightly the boundary conditions obtained in this section to account for the curva- 
ture and to employ less direct arguments to demonstrate rigorously their consistency with 
the bulk gauge condition. 

A. Reflection of electromagnetic waves off a planar perfectly conducting bound- 
ary in Lorentz gauge 

We consider classical electrodynamics in Lorentz gauge in a semi-inflnite flat spacetime 
bounded by a planar accelerating perfectly conducting wall. For concreteness we shall 
assume (3-1-1) dimensions, even though the generalization to other dimensions is straight- 
forward. Let z = Zwit) be the trajectory of the moving boundary, assumed timelike but 
otherwise arbitrary. In the semi-inflnite region z > Zu,{t), with t, x, y arbitrary, the wave 
equation 

A^,u, = (3) 

describes the forward time evolution. There remains, however, some residual gauge free- 
dom, which may be completely characterized by the scalar solutions to the wave equation 

A,.. = 0. (4) 
5 




FIG. 1: Propagation of the residual gauge ambiguity. Specifying initial data on the 
Cauchy surface S, taken here to coincide with t = 0, z > 0, completely specifies the gauge 
in the diagonally hashed region. However, in the horizontally hashed region, additional gauge 
information is required to specify how the longitudinal pure gauge modes scatter or propagate in 
from the timelike boundary at z = 0. 

It follows that the gauge transformation 

A,^A'^ = A,- A,^ (5) 

preserves the Lorentz gauge condition 

A,,, = (6) 

and describes the same classical field configuration. The residual gauge freedom is quite 
limited. Once initial boundary data have been specified (i.e., on a past Cauchy surface), no 
ambiguity remains as to the evolution of the vector potential forward in time. Once A^ and 
its normal derivative have been specified on a Cauchy surface, the future propagation is 
completely unambiguous, because the initial data for A may be propagated forward in time 
from the same Cauchy surface. For a spatially semi-infinite region, however, pure gauge 
modes may propagate in from the timelike boundary, as indicated in Fig. ^ In order to 
fix the time evolution completely, some sort of gauge condition on this timelike boundary 
is required to specify what longitudinal modes, if any, propagate in and how longitudinal 
modes are reflected off the boundary. 

We first consider a boundary moving at a constant velocity, which may be taken to 
be situated at z = 0. E and B vanish inside the conducting wall. (We assume that 
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the conductor has no magnetic flux frozen in.) It follows that Ey and vanish at the 
boundary, although E_l and By need not vanish there, because these may be generated 
by surface charge and current density, respectively, on the boundary. We flx the gauge 
inside the conductor by setting A = {Aq,A) = there. It follows that Ay and Aq — or, in 
other words, all the components parallel to the boundary — must vanish. A jump in A_i_, 
however, is not excluded. Such a jump can be generated from a gauge transformation of 
the E = B = conflguration, for example one where A vanishes inside the conductor and 
rises linearly with distance from it on the exterior. 

In the bulk z > 0, for flxed wavenumber k = (ky, k±), there are three modes consistent 
with the Lorentz gauge condition in eqn. two physical modes 

e|^-^)(k) exp[+ik-x], ^12){^) exp[+ik-x] (7) 
and a longitudinal mode 

ef^)(k) exp[+zk-x], (8) 

where 

VI -(k -6^2 Vi_(k.ey)2 
and k = k/|k|, = (1, 0, 0, 0). A fourth polarization 

e(F) = ^(et-k), (10) 

which we shall call the forbidden polarization, is excluded by the Lorentz gauge condition. 

We have already specifled boundary conditions for the components of the 4- vector poten- 
tial parallel to the conducting boundary, in the present special case At, A^, and Ay. These 
components satisfy Dirichlet boundary conditions. It remains to determine the boundary 
condition for the normal component A^, A^ for the special case here. It is apparent that for 
incoming longitudinal modes to scatter off the conductor into outgoing longitudinal modes, 
without creating any of the forbidden fourth polarization, a Neumann boundary condition 
is required, in other words, Az^z = 0. More generally, 

t-A = 0, Vn(n-A) = 0, (11) 

where t is an arbitrary vector tangent to the conducting surface and h is the normal vector. 
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This nice separation into physical and pure gauge (longitudinal) modes is not preserved 
under Lorentz boosts. While a pure gauge mode always transforms into another pure 
gauge mode, a physical mode transforms into a physical mode with some, in general non- 
vanishing, admixture of the longitudinal (pure gauge) mode. It is not possible for observers 
moving relative to each other to agree on a common notion for the absence of a longitudinal 
component. 

For the case of a stationary conducting boundary, it is possible to remove this ambi- 
guity by further fixing the gauge to Coloumb gauge, where Aq — 0, which is equivalent 
to postulating the absence of longitudinal photons for the preferred stationary observers. 
However, when the velocity of the boundary is not constant, physical photons reflected from 
the boundary will, upon being diffracted back, acquire a longitudinal component from the 
point of view of an observer at rest with respect to the boundary at a later instant. Conse- 
quently, it is necessary to admit longitudinal modes when considering curved or accelerating 
boundaries. 



B. Reflection off a planar perfect gravitational mirror 

We now consider linearized gravity for a metric perturbation h^^ of Minkoswki space 
having the background metric ri^i,. Under a gauge transformation generated by the displace- 
ment field ^'^(x), the metric perturbation /i^,^, defined so that g^^,^ = 77^,^ -|- h^,^, transforms 
as 

hfj,i/ > /i^jy = h^i/ — ^n^i/ — (12) 

It is convenient to define the trace-reversed metric perturbation 

hf^u = h^^ - ^Vni^V'^'hpa. (13) 

[In our tensor calculus, valid only to first order, indices are raised and lowered with the 
unperturbed metric rjixv] We impose the Lorentz gauge condition 

Ku,. = 0. (14) 

It follows that a displacement field satisfying 

□ Cm = (15) 



8 



preserves the gauge condition in eqn. (jl4j) . and that the equation of motion for the metric 
perturbation becomes 

□ V = -(16vrG) T^,, (16) 

which for a vacuum stress-energy tensor becomes 

□ Ku = 0. (17) 

If we consider plane waves with the spatial-temporal dependence 

exp[2k ■ x] = exp[ik ■ x — iut], (18) 

we find that there are a total of ten possible polarizations for h^,^: two physical polarizations, 
four longitudinal (pure gauge) polarizations, and four forbidden polarizations, excluded by 
the Lorentz gauge condition. Explicitly, these are as follows. The two physical polarizations 
are 

h+) = ^(e(i)®e(i) -e(2)®e(2)), 

kx) = ^(e(i) ® 6(2) +6(2) ® e(i)). (19) 

The four longitudinal (pure gauge) polarizations are 

e(9F) = -^(e(i) 6(1) + e(2) ® 6(2)), 
e(gL) = %) ® e(i), 

= ^(e(i) ® %) + e(L) ®e(i)), 

^(92) = -^(e(2) ® e(i) + e(i) ®e(2)), (20) 
and the four forbidden polarizations are 

^ifF) = e(F) ® e(^F), 

kfL) = |, 

%i) = ®e(i.) + e(^) ®e(i)), 

km = ^(e(2) ® e(i.) + e(^) (g)e(2)). (21) 

Here the double hats denote unit vectors for a second-rank tensor. 

As a boundary condition, let us first consider an ideal gravitational mirror in the form 
of a brane with vanishing stress-energy density localized on the brane and a Z2 orbifold 
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symmetry across the brane. We take the geometry of the unperturbed brane to be flat. 
This ideahzed invisible brane is the short distance hmit of any sort of actual brane, because 
for any actual nonsingular brane, there is a finite length scale characterizing its vacuum 
curvature and that inducing by the matter on it. 

In the usual Cartesian coordinates, we obtain the following boundary conditions for a 
normally incident gravitational wave. Let ti and t2 be arbitrary vectors tangent to the 
brane. One has the boundary condition 



as derived from the Israel matching conditions, where d/dn denotes the inward normal 
derivative. We deduce the boundary conditions for the remaining components by consid- 
ering the scattering of incoming longitudinal (pure gauge) graviton modes. It is necessary 
that these scatter exclusively into outgoing longitudinal modes. It is evident that the 
boundary conditions 



for the other components satisfy this requirement, for normally as well as for obliquely 
incident waves. 

III. LORENTZ GAUGE IN A CURVED BACKGROUND SPACETIME 

In a curved spacetime, under the gauge transformation generated by an infinitesimal 
(linearized) displacement field the linearized metric perturbation transforms according 
to 




(22) 




(23) 



hab 



h'ah — hab — ^i9ab " " Ca;6 " 6;a- 



(24) 



The trace modified metric perturbation consequently 




(25) 



transforms as 



h'ab — hab - Ca;6 " 6;a + Qab C;c- 



(26) 



10 



For gauge transformations that preserve the Lorentz gauge condition eqn. Shab^b = 0, 
which imphes that satisfies 



(27) 



where □ = g^'' VaVf,. The evolution equation for hat is [see Appendix |X1 for a derivation] 



ab 



giThda + g^^Jhb + g!!i%, - g^Thab B!fl + 2i?il = -(IGttG) T^IK (28 



When the background is maximally symmetric, a number of simplifications result. For a 
(d+l)-dimensional maximally symmetric background spacetime 



(0) 



'-bed ^^d^l^ 

SO that eqn. (PTj) becomes 



i/ac gbd gad Ubc I 5 



R. 



(0) 



(0) 



(29) 



□ 



d+1 



(30) 



and eqn. (^Hjl becomes 



rf+ 1 



-(rf- l)hab + gfb h. 



-{IQtxG) 



(1) 



(31) 



For AfiS*^ presented using the Randall-Sundrum line element 

1)2 



dz^ - dr + dyi' 



(32) 



we find that 



abed 



-77 {gae gbd — gad gbe) , Rab 



-p9ab, R 



4-5 



(33) 



and Tab = —^gab where A = —3/(4(7^^), which when substituted in eqn. (jSH) yields 

UKb = 0. 



(34) 



In (4+1) dimensions the graviton field hab has (5 ■ 6)/2 = 15 independent components 
or polarizations. Under the Lorentz gauge condition five of these are rendered forbidden 
polarizations. The residual gauge freedom may be characterized completely by the solutions 
to the homogeneous vector wave equation eqn. ()27|) . Here ^'^(x) is simply an infinitesimal 
displacement field generating a coordinate reparameterization. This residual gauge freedom 
corresponds to five longitudinal (pure gauge) graviton polarizations, leaving a total of five 
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"physical" polarizations. Again, observers moving at a constant velocity with respect to 
each other cannot agree on a common notation of a purely physical mode without any 
admixture of longitudinal polarizations. 

Before proceeding to a further analysis of the residual gauge freedom in Lorentz gauge, 
let us first, for comparison, review the counting of allowed polarization states and residual 
gauge freedom in the Randall- Sundrum gauge (see Ref. , where h^^ = /is^ = /i^^ = 

h^u,u = 0. There are five remaining polarizations, all of which are physical. Fixing the gauge 
completely in this way in general displaces the brane from its non-perturbed position at 
z = 1. Consequently, in order to account for all the physical degrees of freedom, it is also 
necessary to include an additional scalar field, which we shall denote ^brane±{^), localized 
on the brane, indicating a normal displacement of the brane with respect to the bulk. [In 
our notation, the vector x denotes a position on the brane whereas x denotes a position in 
the bulk.] In Lorentz gauge this last degree of freedom is absent because the corresponding 
physical degree of freedom is instead described by one of the five longitudinal graviton 
polarizations. If we dispense with the normal displacement of the brane with respect to 
the bulk geometry, then the brane is fixed to the surface z = -2(6) (t), with t, xi, X2 and X3 
arbitrary, no matter what the perturbation is. 

Consider now a gauge transformation generated by the displacement field ^"(x)- For this 
to be pure gauge, we must displace the brane about its unperturbed position z = -2(6) (t) 
according to the normal displacement field ^6rane±(x) = na^''(x). However, if instead we fail 
to compensate for the change in position of the brane by setting ^6rane±(x) equal to zero, 
what would be a pure gauge transformation corresponds to a displacement of the brane 
trajectory relative to the bulk geometry. The remaining four longitudinal polarizations 
allowed in Lorentz gauge correspond to reparameterizations of the brane. 

The Lorentz gauge condition in the bulk does not in any way fix the gauge on the brane. 
A gauge transformation corresponding to a reparameterization of the brane is effected by a 
flux of incoming longitudinal gravitons emanating from the Cauchy surface and reflecting off 
the brane. More specifically, let the field (x) indicate an infinitesimal reparameterization 
of the brane. The longitudinal gravitons implementing the desired gauge transformation are 
generated by choosing ^(x) on the initial Cauchy surface such that ^||(x) on the boundary 
is reproduced and ^±(x) vanishes there. Because we do not include an additional degree of 
freedom corresponding to normal displacements of the brane, from the view of an observer 
on the brane, the longitudinal modes in the bulk giving vanishing 'C||(x) but non- vanishing 
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'C±(2£) on the brane are not "pure gauge". These longitudinal modes take the place of the 
field localized to the brane ^brane±{l^) present in the Randall- Sundrum gauge. 

IV. COMPATIBILITY OF SOURCES ON THE BOUNDARY WITH THE 
GAUGE CONDITION 

In this section we investigate the compatibility of boundary conditions on the timelike 
boundary (i.e., the brane) with the Lorentz gauge condition. First, as a sort of warm-up 
exercise, we consider a fiat bulk with a planar boundary, which acts as a perfect gravita- 
tional mirror. Then we consider the modifications necessary to take into account extrinsic 
curvature (i.e., nonvanishing stress-energy on the brane at zeroth order) and the spacetime 
curvature of the bulk at zeroth order for perfect AdS. 

A. Special case — the planar perfect gravitational mirror embedded in a 
Minkowski bulk revisited 

We first consider the following boundary conditions for a linearized source on the gravi- 
tational mirror boundary (here, in this special case, the bulk is Minkowski space and so is 
the unperturbed brane) 

hAB,N = —i^%,)Tab-, (35) 
hAN = 0, (36) 
hNN,N = 0. (37) 

Here the indices A,B,... indicate directions tangential to the brane, the normal direc- 
tion, and a,b, . . . all five directions. Tab denotes here the linearized stress-energy pertur- 
bation on the brane. Eqn. ()35|1 is the Israel matching condition and eqns. ()3fi|l and (jTTj) 
are supplementary gauge conditions on the boundary. We show that the Lorentz gauge 
condition is subsequently satisfied if and only if it is satisfied on the initial Cauchy surface 
and Tab,b = as well. 

To this end, we consider the vector field 

Va = hab,b, (38) 

assumed to vanish and have vanishing normal derivative on the initial Cauchy surface. If 
each component of this 4-vector field can be shown to satisfy either Dirichlet or Neumann 
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homogeneous boundary conditions on the timehke boundary (i.e., the brane), it follows that 
the Lorentz gauge condition holds everywhere. Physically, this means that wave packets 
obeying the gauge condition reflect into wave packets obeying the gauge condition and that 
any inhomogeneous source on the timelike boundary does not emit any waves violating the 
Lorentz gauge condition. 

On the boundary, the relation 

Vn = hNN,N + hAN,A = (39) 

holds trivially because of the supplementary gauge conditions (jHUj) and (jHTj) . Similarly, the 
normal derivative of the tangential components is given by 

Va,N = hAN,NN + hAB,BN = —hAN,BB ~ I^\,)TaB,B = —l^%yI^AB,B- (40) 

Consequently, no Lorentz gauge condition violating waves emanate from the boundary if 
and only if 

Tab^b = 0, (41) 
in other words, if Tab is a conserved tensor fleld on the brane. 

B. The general case — bulks and branes curved at zeroth order 

In the previous section we were able to show in a few lines, relying on the wave equation 
in the bulk and stress-energy conservation on the brane, that (jHBj) and (jHTj) are the correct 
auxilliary boundary conditions in Lorentz gauge if one assumes that both the brane and 
the bulk are flat at zeroth order. In this section we essentially repeat the derivation of 
the previous section for the case where the boundary and the bulk geometry are not flat, 
including modiflcations necessary to account for the non-vanishing curvature. When the 
bulk is curved, the condition Jian = no longer necessarily implies that hAN;A = on 
the boundary. Consequently, it is necessary to modify eqn. (P7|) to contain non-derivative 
contributions of the metric perturbation. 

To treat the curved case, we flrst flx a point x on the boundary E. Then we coordinatize 
a neighborhood of x on S using the coordinates xa chosen to be Fermi normal coordinates — 
that is, using exponential map from the tangent space of x on S to S. By choosing Fermi 
normal coordinates we set all the Christoffel symbols of the metric connection on S to zero 
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at X, although the first and higher partial derivatives of the Christoffel symbols do not in 
general vanish there. These coordinates may be continued off the boundary by means of 
the Gaussian normal prescription — that is, the coordinate n in the N direction vanishes 
on E and the coordinates on E are continued along initially normal geodesies, n indicating 
physical distance along these geodesies from E. In this way, at x we obtain the following 
connection coefficients for the unit vectors ca — d/dxA and ejv = d/dn : 

VatCat = 0, 

V^eAT = +Kab es, 
VatBa = +Kab gb, 

VaCb = -Kab eN, (42) 

where the tensor Kab is the extrinsic curvature and V denotes the bulk covariant derivative 
induced by the bulk metric. These equations may be explained as follows. The first equation 
simply defines the continuation of N off the brane. The second equation consists merely a 
restatement of the definition of the extrinsic curvature. The third equation follows from the 
vanishing of the Lie derivative of along the brane (i.e., [^, iV] = 0). Finally, the fourth 
equation follows because is orthogonal to E. Away from x, the last equation would be 
modified to 

VaCb = i^AB ec - Kab cat, (43) 

where F^g denotes the connection coefficients of the covariant derivative on E with respect 
to the metric induced on E by the bulk metric. For future reference, we also give the 
following second covariant derivatives along the brane 

VdVc Cat = +Kca;D ^a - Kca Kda cat, 

VdVc cb = (f cs;£) - KcB Kda) ga - Kcb-d gn- (44) 



1. The electromagnetic case 



For the electromagnetic case we calculate on the boundary at x 



— ■A^-N' + 



[Ag) CB + Ai^) 



= An-n + Va 

—^ + KaaA^^K 

OXb 



■ eA 



(45) 
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Here the superscript (c) denotes that the component in question is to be regarded as a scalar 
quantity (i.e., differentiated using ordinary partial rather than covariant differentiation). 
From our boundary condition Ab = (i.e., t- A = 0), it follows that the second term 
vanishes. However, to ensure the reflections and sources from S do not violate the Lorentz 
gauge condition, it is necessary to modify the supplementary boundary condition to 

An;N + Kaa An = 0, (46) 

so that eqn. P5|l is set to zero. 



2. Linearized gravitational perturbations 

For linearized gravity on a curved boundary only the following boundary condition stays 
unaltered 

hAN = 0. (47) 

This is because adding any derivative term would spoil the short distance limit. Conse- 
quently, we retain the boundary condition in eqn. (jTTj) and next find the modifications to 
eqn. (|?fj) required for the curved case. We find that 

Vn = hNN;N + hNA;A 

+Va h^NN i^N ® Cn) + fl^ND i^N ® 6/)) 

+hcN (ec ® Cat) + Ji'^cd (ec ® cd) ■ {cn ® ca) 
= hNN;N + hjyN Kaa — hsA Kab- (48) 

The above implies that the boundary condition valid for the case with a fiat brane and a 
fiat bulk, hNN,N = [eqn. (jHTf) ]. must be modified to 

hNN;N + Kaa h^N - Kab hAB = (49) 

for the general case. At this point that the boundary condition eqn. (jTfjl remains unaltered 
is just an inspired guess, vindicated in the calculations that follow. 

We now compute the normal derivative of the tangential component of the gauge con- 
dition, 

Va-N = hAN\NN + hAB:BN 
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— hAN;BB + hAB;NB 
+ {hAN;NN + hAN;BB) 
+ {hAB;BN — hAB;NB)- 



(50) 



The second line of the second expression can be simphfied by using the wave equation in 
eqn. and the Gauss-Codacci relations in eqn. (jBljl 



hAN;NN + hAN;BB — D ^AAf — —'2'RnDAC hcD — "2 (^K^DA;C ~ ^DC;a} hcD, (51) 

and the third line simplifies using 

hAB;BN — hAB;NB = RbCNB hAC + RaCNB hcB = {KbA;C ~ KbC]A) hcB, (52) 



so that 

Va;N = —hAN;BB + hAB;NB — (-^b1;C ~ ^BC;a} ^cB- 

We now evaluate Tian-^bb- From the definition of the Laplacian operator 

= (iV^V^) (iV-V,) - {N^^V^N-) V^ + Y. [(A^V^) (A-V,) - (A^V^A-) VJ 

A 

= Vat Vat 

A 

it follows that 



(53) 



(54) 



hAN: 



BB 



VbVb + K'^BB^n] hAN. 



(55) 



To simplify Han-^bb exploiting the fact that Han = on the brane, we compute 



^D^chAN 



VnV 



+h!'EN i^E (S) Cat) + h'-^^p {ce ® ep) ■ (ca ® Cat) 
Voh'^NN (VcCAf) ■ eA + ^ch^^N (VdCat) ■ Ca + /ij^AT (VdVcCat) ■ ca 

+^AF (VdVcCf) ■ eN + ^Dh^AF c^f) ' + ^ ch% (VdCf) ■ Cat 

J^AC '>'NN;D + J^AD i^NNiC + J^AC;D ^^NN 
-K^FC-D hAF - K^F hAF;D " -K'df /^AF 



';C- 



(56) 



Contracting C and D, we obtain 



V B^ BhAN — '^K^AB hNN;B + K^AB;B ^NN — K^C^;B ^AC ~ '^K^BC ^AC^B- 



(57) 



17 



Consequently, 

Va-N = hAB;NB — hAN;BB — (-^iM;C ~ ^BC;a} ^CB 

= hAB;NB — '^K^AB hNN;B + '^K^BC hAB;C — K^'^^ hAN;N 

-K^Z,B hNN + kPc;C hAB - [K'-iXc - K^ca) hcB. (58) 

It remains to be shown using the Israel matching condition and stress-energy conserva- 
tion on the brane that the above quantity vanishes, so that the corresponding boundary 
condition on the brane is homogeneous. 

V. STRESS-ENERGY CONSERVATION ON THE BRANE AT FIRST ORDER 

In this section we consider how to express stress-energy conservation for a singular 
distribution of co-dimension one of stress-energy localized on the boundary brane about 
which a Z2 symmetry has been imposed. The simplest way to proceed, which does not 
require any additional principles for how to treat the singular distribution of boundary 
matter, is to consider in the limit 5 — > -|- a symmetric distribution of non-singular Z2 
symmetric matter of a finite small thickness 25. The bulk spacetime is reflected about the 
hypersurface at the center of this boundary stress-energy by means of the Z2 symmetry, as 
indicated in Fig. |21 The four-dimensional projected stress-energy tensor on the brane Tab 
is obtained by integrating over the coordinate normal to the brane, denoted by A^, whose 
units correspond to physical distance. One has 

Tab{^) = r'dN Tab{^,N), (59) 
J —S 

where Tab is the five- dimensional stress-energy tensor. For finite 6 this projection procedure 
is not entirely satisfactory because of ambiguities in parallel transport. However as 5 — > -|- 
these difficulties disappear. 

Using ordinary derivatives, we may write 

dN {Tab,b + Tan,n) = Tab,b + Tan{N = +S) - Tan{N = -S) 

= TAB,B + 2TANiN = +5). (60) 

Here we may choose, for the purpose of this calculation, Gaussian normal coordinates, so 
that = corresponds to the hypersurface S of Z2 symmetry, and is the signed normal 
physical distance from E. In the neighborhood of a given point x G S, we may choose the 
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FIG. 2: Schematic diagram of the Z2 symmetric distribution of brane stress-energy. 

A singular distribution of stress-energy on the brane is most easily regarded as the limit of a 
non-singular distribution of stress-energy of support confined to within a shell of thickness 25 in 
the limit 5 — > -|- . We derive the boundary conditions applicable for an infinitely thin singular 
brane by taking the 6 0+ limit of a nonsingular thick brane of thickness 6. The (dashed) line 
into the middle represents the Z2 reflection orbifold symmetry. The projected stress-energy Tab 
(always tangent to the brane) is distributed over the indicated shaded region of thickness 26. 

transverse coordinates, labelled A,B,..., so that the connection coefficients vanish at x. 
If 6 is sufficiently small, the difference between (Tab,b + Tan,n) and (Tab;B + Tan;n) is 
negligibly small, and it follows that in the 6 — > 0+ limit 

Tabib + 2Tan{N = +5) = 0, (61) 

where | indicates the covariant derivative on S, or equivalently (in the 5 0+ limit) on 
the displaced surface (S + 5) just outside the singular distribution, on which one imposes 
the Israel matching condition. Here —Tan represents the five- dimensional flux of energy- 
momentum incident on the brane from the bulk (the direction is taken to be outward), 
or with a plus sign it would represent energy-momentum flowing from the brane into the 
bulk. 

We note that, alternatively, we could have derived eqn. (pT|) directly from the Israel 
matching condition and the Gauss-Codazzi relations in eqn. ()B1|) . 

Because we have assumed an empty bulk (except for the negative cosmological constant), 
with Tab = —^''^^dab outside of the brane, the flux contribution 2Tan in eqn. (16 Ij) vanishes. 
This follows because the direction A is normal to A^. Extra degrees of freedom (e.g., bulk 
scalars) in the bulk would in general spoil the vanishing of Tan- 

In the discussion above, the (projected) brane stress-energy Tab has been shown to 
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be a tensor field parallel to the brane in both indices and conserved with respect to the 
four-dimensional metric induced on the brane. For the components parallel to the brane, 

Tab\b = Tab;B = 0, (62) 

where | and ; denote covariant derivation with respect to the four- dimensional metric in- 
duced on the brane and the full five- dimensional metric, respectively.^ The discussion so 
far in the present section applies equally well with and without perturbations. 

We now consider perturbations to hnear order, splitting Tab — T'ab + '^ab- zeroth 
order, stress-energy conservation on the brane is expressed as T^^^^ — 0, where | is the 
connection with respect to the unperturbed metric induced on the brane. 

Stress-energy conservation to linear order is obtained by extracting the first-order terms 

(,fof - /.-^) vT^'^ (rig + rw) = 0. (63) 

Here the superscript [^f^^^ -|- h] indicates that the covariant derivation is with respect to the 
perturbed metric. Keeping only the first order terms of the above equation gives 

_ I.BC ^(0) ^[h] ^(0) _ „ .X 

^AB\B ^AC\B^^B -^AB—^J 

where | denotes the connection of the unperturbed metric induced on the brane. Here the 
operator v|^'' denotes the change in the metric connection due to the perturbation of the 
metric from g^^^ to (g^^^ + h), accurate to linear order in h. We may express this operator 
as 

Vjles = ^AB ec, (65) 

where the pseudo-Christoffel symbols above constitute a third-rank tensor (and not a 
pseudo-tensor) with respect to the unperturbed metric, and can be expressed as 

r''*'cAB — 2 hAC\B + hBC\A — hAB\C ■ (66) 

At a given point this relation can be demonstrated by choosing Fermi normal coordinates 
there with respect to the unperturbed metric. Precisely at this point and in these coordi- 
nates, g'AB,c — pl^'"'] '--^ AB — 0. It follows from the generally valid expression 

r^AB = ^5''^^ 9ad,b + 9bd,a - 9ab,d , (67) 



^ Note that, for example, with the five-dimensional metric connection Tan;B does not necessarily vanish, 
even though Tan is zero; however, in the equations these terms are not included. 
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that to linear order 



AB - 7j5'(0) 



hAD,B + hBD,A — hAB,D 



(6^ 



Because the unperturbed Christoffel symbols all vanish, these ordinary derivatives may be 
replaced by covariant derivatives, rendering equation eqn. (j68p equivalent to eqn. . It 
follows that 



^B ^AB - -i ' ' BA -IcB - ^ BB i AC 

1 ' 

"2 



= —-j{hBC\A^hcA\B — hBA\c)^CB 

-^(^h^dB + hcB\B — ^-B-B|c)^C 
1 _ , . 



9-r(0) . ^(0) , ^(0) , 

^IaC 'T'BCIB — -IaC '>'BB\C + -^CB '>'BC\A 

The first-order terms of the stress-energy conservation equation give 

[yB-lAB] - ^AB\B ~ '^BC ^ AC\B + ^ B ^ AB 

1 



- ~ '^BC ^AC\B ~ ^AC I^BC\B " -^BC riBC\A + --i AC 



0. 



(69) 



(70) 



VI. THE LINEARIZED ISRAEL MATCHING CONDITION AND ITS DIVER- 
GENCE 

In this section we show how to derive the form of the linearized Israel matching condition 
for a Z2 symmetric stress-energy distribution of co-dimension one in a (4+1) -dimensional 
bulk spacetime. The exact Israel matching condition 

,^2 



Kab 



"(5) 



'^ab — -Qab'^cc 



gives K = H — fp-'?', which may be used to re-write eqn. ()7H1 as 



(71) 



Kab - Qab K 
giving us at first order the equation 



"(5) 



%iB, 



K^l - 9% {K'^l - hen) - hAB 



K 



(5) ^(1) 
TT ■'■AB- 



(72) 



(73) 



The first order perturbation in the extrinsic curvature is [see AppendixOfor a derivation] 

1 1 

= 2 \-hAB;N — hAN;B — hBN;A] + 2^^^ ^nn- (74) 



K 
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Using the supplementary boundary condition h^N = 0, we may simplify the terms hAN;B + 
hBN;A- Expanding 



'l'AN;B — V_B 



h% Cat ■ (VbCd) + h^§j^ ca ■ (VbCat 



(75) 



we may now re-write 



K 



(1) _ 1 



AB 



hAB;N + K^BC ^AC + Kac hsc 



-K^'l h 



AB "-NN, 



(76) 



from which we find that 



^CC — 2 CC;N + J^CD '''CD — 2 C-C* "-A^Af- 



(77) 



The next step is to re-express eqn. fj76|) in terms of the trace-modified metric perturbation 
hab rather than hat- From 

1 



it follows that 



hAB — hAB — ^dAB {h-CC + ^A^Af) 



hAB-N — hAB;N — TtS'aB {hcC;N + hNN;N) 



We use the boundary condition [from eqn. PH)) ] 

hNN;N + Kaa hNN — K^AB ^AB = 



(7J 



(79) 



^0) 



and 



hNN;N — hNN;N — -Qnn {hcC\N + hMN;N) 



{hNN;N — hcC]N) (81) 



to obtain 



hNN;N — hcC;N + '^hNN;N — hcC;N — "^Kaa + '^Kab hAB, 



^2) 



which gives 



hAB;N — hAB;N — QaB hcC;N + 5'aB 



,(0) 



^CC ''NN — J^CD '''CD 



Using eqns. (fTSj) and (jTHl) and the relations 



- 1 - 2 - 

hab = hab - ^9abh, h = --h, 



^3) 



^4) 
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we expand the left-hand side of the perturbed Israel matching condition in eqn. (j73j) as 



K'^l - afs {kSI - K§1 hen) - K^'' h^s 
\ {hAB;N + k'^I hsc + Kfc hAc) - \k'^b h 



NN 



(0) 

-9ab 



^hcC;N + Kgl hcD - h 



NN 



K^D hcD - K^'^ h^s 

2 {hAB;N - 9aB hcC;N) + ^ {K 

\{Kfs-9fBK^'')hNN-K^'h 



AC '''BC + ^BC '''AC 



2 
1 r 



hAB;N + 9% {K^cI hcD-K^'^ h 



AB 

NN 



+ 2^ic [hBC - ^9Z h) + -KZ [hAc - ^9Z h 



.(0) 



,(0) 



1 



h 



NN 



11 1 

-hAB;N + 2^ AC ^BC + 2^BC ^AC 

h) . 

3 / 3 

'\kP^ hcD - \k^'^ hNN + Ik('^ h 

\hAB;N + Ik^ac hBC + IkPc hAc - K^'^ h 

,1(0) 7^(0) T 1^(0) T 1^(0) 1,1(0) r^{0)T 

+ 2^AB ^CD "-CD - 2 AB "'NN " AB + q9aB J^' 'n. 



'(0) 



1 r 
2 



kTb-9%K^'' 



2 

hNN 



K^'^ hAB 

1 ^.(0) 
^AB 



, (0) 

+9ab 



i-AB 



^5) 



We re- write eqn. (I73|l . obtaining for the linearized Israel matching condition 



-iT-ABiN + 2 '^BC 
' ,(0) 7^(0) 



+ Ik^I h^c - K^'^ h^B 



+ 2^AB hcD - 2^AB ^NN 



"(5) ^(1) 
■^AB- 



6 6"^'' 2 

We now proceed to take the divergence of eqn. (j86p considered as a tensor field 
brane using the covariant derivation V generated by the induced metric on the 
Noting that 



(86) 

on the 
brane. 



N 



hAB;N\B 
hAB;NB - 



1^(0) I 7^(0) T 7^(0) , 

J^AR n-NB-.N — J^RR riAN:N + J^nR ll 



CB "'AB;C, 



(87) 



we find the divergence of the left-hand side of eqn. (|H^ is 



B 
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1 



AB:NB 



AB;B — Ti^AB 'T'NNiB 



2AB iT'NBiN — l^^BB I^AN-^N + l^^CB I^ABfi 
+ ]^Kfc hAC;B + \Kfl hBC;B + \Kfc hBC;A ' K^'^^ h 
-Ik?) + Ik^'^ kA 

+ ^KbC;B hAC + ^^ic;B ^30 + ^^BC;A ^BC " K^B ^AB " ^K^AB;B W 



— K 



(0) 

AB:B 



h + ^k5^ h . 



6 6 

We next compute the divergence of the right-hand side of eqn. ()86|) . using the result from 
the stress-energy conservation at first subleading order [eqn. (|7nj) ] 



r: 



(1) 



AB\B 



T, 



(0) 



AC\B "'BC 



V R i / 



AC|B "'BC' ^AC /^BCIB + -^sc f^BClA - --Iac '^BB\C- 



^9) 



Because of the boundary condition Kan 
matching condition 



0, hAB\c = hAB;C- The zeroth-order Israel 



K 



(5) ^(0) 
^AB 



(90) 



allows one to express eqn. (IHUjl as 

,2 



15^^(1) 

2 -'A-BIb 



if 



(0) 

AC]B 



(0) 7^(0) 



B 



iBC 



+Kfl hBC;B - K^''^ hBA;B + \Kfc hBC;A ' i;K^AC hBB;C- (91) 



1 

2' 



Using the definition of the trace-modified metric perturbation [eqn. (j84j) ]. we find that 

,2 



K 



i51^(l) 

2 ■'abib 



^AC;B - 9 AC J^;B [ I^BC - ^9bC " 



h 



(0) j^(o)- 

/iBA;B 



AC "-BCiB 
+ M°<^/^;C+ V) ^A, 



-A^C* /^BCiA - 2 AC '^BB;C 



3 . g 

and noting that h = Jibb + we can re-write eqn.(|U^ as 



(92) 



K 



2 ■'abib 



kTc;B - 9Z Kr^' ( hBC - ^9Z h 



.(0) T^m 



+KZ hBC 



1 

6' 



-ii^i") h,c + /^;A. 



1 

6' 



/iBA;B + ^i^g'i hBCA + ^^AC ^^^;C 



(93) 
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Subtracting eqn. (j93|) from eqn. (j88p we find that vanishing of the divergence of the Israel 
matching condition and stress-energy conservation on the brane imply 



+ -hAB;NB 

-2^AB 'T-NB;N " 2 BB '^AN;N " 2 AC '^BC;B 
'^^BC ^AC;B — -^AB h^NiB 

+ 2^BC;B riAC - 2'^AC;B '^BC + -j^i^BC.A I^BC - -^J^ AB;B I^NN 

+\KfB;B h - \k^'2 h = 0. (94) 

The zeroth order Israel matching condition [eqn. (jHOI)] and stress-energy conservation at 
zeroth order T^^^^ = yield 

K%.,B - = 0, (95) 
which may be used to rewrite eqn. ()94p as follows 



2^AB;NB 

1 



AO) h J- T ^ ^(0) 7 



2 2 ^iY,iv 2 

I lAo) I lA^) u 

-i-n^BC '>'AC;B — J^AB '''NN;B 

+ \Kfc;B hAC - lKfc;B hBC + ^i^^^^A W - ^K^^-^B W = 0. (96) 

We note that eqn. ()9(ij] resulted from assuming the Israel matching condition and stress- 
energy conservation on the brane. Our object was to show that Va-n, given in eqn. 
vanishes. We note that eqn. equals precisely one half the right-hand side of (jSHl)- 
Hence it follows that Va-n = 0. 

VII. DISCUSSION 

We summarize our main results as follows. Throughout the paper we have assumed 
Lorentz gauge in the bulk for the linearized metric perturbations — that is 

hab;b = 0, (97) 

where in terms of the linear metric perturbation hab 

hab = Kb-\gf^ {g^'^^'Ka). (98) 
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In the bulk the Lorentz gauge condition is consistent as long as the bulk stress-energy is 
conserved. We have demonstrated that the auxiliary boundary conditions on the brane 

hAN = (99) 

and 

hNN,N + Kaa hNN — Kab Tiab = (100) 

together with the Israel matching condition sourced by a conserved stress-energy on the 
brane Tab ensure that the waves reflecting from the boundary or emanating from sources 
there respect the bulk Lorentz gauge condition. 

The bulk wave equation in (4 + 1) dimensions [given in eqn. ()28|1 . which simplifles to eqn. 
flH4|l for an AdS^ bulk] propagates 15 components. On the boundary the Israel matching 
condition provides 10 of the 15 required boundary conditions. Eqns. ()99|) and ()100|) provide 
the remaining 4+1 boundary conditions, respectively, required to render reflection off the 
boundary unique and consistent. 

The Lorentz gauge condition in the bulk [eqn. (j97j) ] does not flx the gauge on the brane. 
Rather, once a choice of gauge on the brane has been chosen, it provides a prescription 
for continuing this gauge choice into the bulk. In (4+1) dimensions the Lorentz gauge 
condition allows exactly flve longitudinal (pure gauge) modes in the bulk. Four of these 
correspond to reparameterization of the brane. The one remaining mode corresponds 
to normal displacements of the brane, which become physical if the brane is flxed to its 
unperturbed position with respect to the unperturbed background spacetime. 
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APPENDIX A: BULK METRIC PERTURBATION EVOLUTION EQUATIONS 
IN LORENTZ GAUGE 

In this appendix we derive the evolution equation for the metric perturbation hah and 



ab 



hOabR 



for the violation of the gauge. We assume the Einstein equation Gab 
{87iG)Tab and the Lorentz gauge condition hab-b = 0. We also demonstrate the consistency of 
the Lorentz gauge condition in the bulk assuming that the bulk stress-energy is conserved. 
For the linearized perturbations 



G 



(1) 



R 



(1) 

ab 



gj^R^'^ + habR^'^ = {8nG)T^ 



,.(1) 

- ab 5 



(Al) 



where i?*-^^ = g^^^ R^^J — h^'^ RfJ- The Riemann tensor is given by 



R" 



bed 



ra "pa I "pa "pe "pa p^ 

bd.c — J- bed + J- ec J- bd — ^ ed ^ be 



(A2) 
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where Vbc = \g'"^ [9dc,b + gbd,c - 9bc,d] can be decomposed as T^^ = rt^'*°'l^^ + r''']"^, 



where T^' 



|fl'(o) [hdc:b + hbd-c — hbc-d] ■ The hnear perturbation of the Riemann tensor 



be 



is R^^^lcd = ^^'^Kd-c ~ '^^'^Kc-d^ for the Ricci tensor 



acb 



- \-h-nh - h. 



ab l^ab;cc ~l~ hca\bc ~\~ ^cfe;ac 



(A3) 



From the definition of hab [eqn. ((23)], hab = hab — {d — 2) ^g^^^h, h = —2{d — 1) ^h. 
Expressing eqn. ()A1|) in terms of hab, we find that 

1 



G 



(1) 

ab 



hab;cc hca;bc ~l~ hcb-ac 9ab ^cd;cd 



-9ab l^cd + riabQcd 



R. 



(0) 
cd 



(SttG) T,«. 



Using hab-cd - hab;dc = Raedc^eb + R^Jdc^ae, and imposing the Lorentz gauge condition gives 



(A4) 



?(o) 



□^afe - gfchda + Qachdb + 9ab Kd " Q^d^ 



.(0): 



R^^ + 2i?il/.,, = -(levrG) Ti^ (A5) 



as the wave equation in Lorentz gauge. 

We next demonstrate the consistency of the Lorentz gauge condition with the above 
evolution equation. To this end, we define the field 



K — hah-b 



(A6) 



quantifying any possible violation of the Lorentz gauge condition and show that Va satisfies 
a homogeneous wave equation — that is, one with no non-vanishing sources, as long as the 
bulk stress-energy is conserved (i.e., Tab^b = 0). The fact that the wave equation for Va is 
homogeneous implies that if Va and its normal time derivative on an initial past Cauchy 
surface vanishes, then Va vanishes everywhere in the future. 
We now take the divergence of eqn. ()A5|) . The identities 



Dab;mn -Dab\nm 
Sabc\mn Sabc;nm 



p r-)C _|_ p n 

^acnm^ b ' ^bcnm^ a ; 

p Q'^ M Q'^_I_R q d 

J^adnm*-^ be ~r -Ttbdnm^a c ' ^cdnm'-^ ab 



(A7) 



give 



'''ab;ccb hab;cbc + Radbchdb;b + Rbdbchad;c + Rcdbchab;d 
= habficc + {Radbchdb + Rbdbchad);c + Radbchdb;b + Rbdbchad;c + Rcdbchab;d 
= \I\hab:b — Radcb;chdb + Rdc:chad — 2Radcbhdb:c + 2Rdchad:c — RdbhabA, 



(A8) 
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which can be used to take the divergence of eqn. ()A5|) 



(nhnK).h - 



9bc '^da;b + Qab '^cd;b\ ^cd + '^^adbc'^dc;b 

gic^h^a + g^lhab + gThcd - gfrhj R^Z + ^R^Xbhdc = -(167rG)T« . (A9) 



Using eqn. ()A8|) . we obtain 



^hab-b — hcd;aR^cd ~ 



/Shdb + gTh. 



cd 



gcd l^ab 
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The zeroth order Ricci tensor can be replaced using 



4? = (SvrG) 



d-1 



gab 



The Bianchi indentity R 



(0) 



d(0) _ p{0) , p{0) , p(0) 
ac[bd;b\ ~ ^bd[ac;b\ ~ ^bdac;b ^Hdcb;a ^bdba;c 



(All) 

gives the relation 



Rbdac:b - RdcL + RdL = 0, which cau be used to re-write 



5(0) , d(0) 
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giving 
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(A13) 



2(rf-i; 

Using the first-order stress-energy conservation equation previously obtained in eqn. (f7n|) 



iVb^afeJ — -iab-M ^ac:b"'bc ^ ac '''bc;b „-^bc ' ''bc;a -r „ -i ac "'bb;c 
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-(1) ^-(0) 



2 - oc ■ ' 2 

■^ab;b ^ac;b"'bc ^ -ibc "-bcia ~^ 2 d — 1 ~ 



(A14) 



we obtain DVa = \3hab-^b = 0, proving the consistency of Lorentz gauge at linear order 
provided that stress-energy is conserved. 



APPENDIX B: DERIVATION OF THE PERTURBED ISRAEL MATCHING 
CONDITION 

We derive the Israel matching conditions for a surface of codimension one having singular 
distribution of stress-energy embedded in a ((i-l-l)-dimensional bulk spacetime. The Gauss- 
Codacci relations js^], 13 

Rabcd = Rabcd + KadKbc — KacKbdi 



32 



RnBCD — KbC;D — KbD;Ci 

Rnbnd = KbcKdc — Kbd,n, 



(Bl) 



decompose the Riemann tensor into a components tangent and normal to the brane. The 
Ricci tensor and scalar are given by 
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Rncnc = KabKab 
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and 



R = RABg^"" + RNNg'''' = R^'""^^ + '^KabKba -K^- 2ir,iv. 



The Einstein equations decompose similarly 



(B3) 



Gab = Gab^ + 2KacKcb - KabK - Kab, 
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SKcdKdc -K'- 2ir,iv 
Ka, 
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(B4) 



where K = g^^a)KAB- 

At the brane itself, where the metric is continuous but its normal derivative may suffer 
a jump, the singular contribution to Gab arises from the terms 



Kab,n + Oab^'K^n, 



(B5) 



which when integrated across the brane give 



Kab - Ks'Tb'^ 
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where we define four- dimensional projected stress-energy tensor 

1-5 



where ^^fd_^_l^ = SnG 



(d+i) — ^"(d+i) 
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Tab = dN Tab, 
M7}\'7?'^ ■ We can trace- modify eqn. ()B6|) to obtain 
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APPENDIX C: DERIVATION OF EXTRINSIC CURVATURE PERTURBA- 
TION 



To compute the perturbation in the extrinsic curvature derived from a perturbation in 
the metric we use the relation 



(CI) 

AB 



This formulation in terms of the Lie derivative circumvents having to consider how covariant 
differentiation is modified by the metric perturbation. It follows that 

1 



%^ + Aiv£^°1^^- (C2) 



For any doubly covariant tensor field a; and contravariant vector field t/, the following 
relation 

(Cjj ^) ^ = W UJab;c + ^ac U^-b + UJcb -a (C3) 

holds regardless of the choice of covariant derivation ; as long as the connection is torsion 
free. This property reflects the fact that the Lie derivative is defined solely in terms of 
the flow induced by C/, and consequently is independent of any metric or affine, torsion 
free structure defined on the manifold. We find it most convenient to set ; to the metric 
connection defined by the unperturbed metric g^^^ ^ in accord with the convention of the 
rest of this paper. To zeroth order one has K^ab — ^a-b- 
We first calculate SN hy requiring that 

(7V« + 5N^) {gf^ + Kb) {n' + 5N') = +1, (7V« + 5N^) {g^^ + h^)x\ = 0, (C4) 

for all directions X\ parallel to the brane. Here X\ is a contravariant unit vector so that 
gl^^ = N'^N^ + Ea^I^a- It follows that 5N^ = =^hNN,SN^ = -/ijv^, which may be 
re-written as 

= 1 [N'KdN'') 7V« - ^("o>6cA^^ (C5) 



so that 



SNa-b = lhNNNa.b+lN'Kd-bN''N, + N'h^N'',bNa-hac;bN'-K 

= ^h^NK^J + lhr,N;bNa + HcnK^^In, - KN;b - KcK'^^I (C6) 
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We now obtain the first order extrinsic curvature perturbation 



-A. 4 o — — 
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